Search for: All records

ABSTRACT Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity $$v_\text{max}$$ and the radius $$r_\text{max}$$ at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for $$v_\text{max}$$ and $$r_\text{max}$$, this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. The technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics. 

Abstract When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature.To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m= 50 cm²/g).When simulating these halos, we also vary the massresolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses.The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10⁶particles per halo), we find that variations of order 10% in collapse time are still present.The results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context. 

ABSTRACT We provide new constraints on the dark matter halo density profile of Milky Way (MW) dwarf spheroidal galaxies (dSphs) using the phase-space distribution function (DF) method. After assessing the systematics of the approach against mock data from the Gaia Challenge project, we apply the DF analysis to the entire kinematic sample of well-measured MW dwarf satellites for the first time. Contrary to previous findings for some of these objects, we find that the DF analysis yields results consistent with the standard Jeans analysis. In particular, in this study we rediscover (i) a large diversity in the inner halo densities of dSphs (bracketed by Draco and Fornax), and (ii) an anticorrelation between inner halo density and pericenter distance of the bright MW satellites. Regardless of the strength of the anticorrelation, we find that the distribution of these satellites in density versus pericenter space is inconsistent with the results of the high-resolution N-body simulations that include a disc potential. Our analysis motivates further studies on the role of internal feedback and dark matter microphysics in these dSphs. 

Abstract We study the evolution of isolated self-interacting dark matter (SIDM) halos that undergo gravothermal collapse and are driven deep into the short-mean-free-path regime.We assume spherical Navarro-Frenk-White (NFW) halos as initial conditions and allow for elastic dark matter self-interactions.We discuss the structure of the halo core deep in the core-collapsed regime and how it depends on the particle physics properties of dark matter, in particular, the velocity dependence of the self-interaction cross section. We find an approximate universality deep in this regime that allows us to connect the evolution in the short- and long-mean-free-path regimes, and approximately map the velocity-dependent self-interaction cross sections to constant ones for the full gravothermal evolution. We provide a semi-analytic prescription based on our numerical results for halo evolution deep in the core-collapsed regime.Our results are essential for estimating the masses of the black holes that are likely to be left in the core of SIDM halos. 

ABSTRACT The gamma-ray Fermi-LAT Galactic Centre excess (GCE) has puzzled scientists for over 15 yr. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology – potentially arising from the annihilation of dark matter (DM) particles – with a boxy morphology – expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VISTA Variables in the Via Lactea survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a galprop-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology for the GCE, regardless of the mask applied to the Galactic plane. 

ABSTRACT We study the evolution of isolated self-interacting dark matter halos using spherically symmetric gravothermal equations allowing for the scattering cross-section to be velocity dependent. We focus our attention on the large class of models where the core is in the long mean free path regime for a substantial time. We find that the temporal evolution exhibits an approximate universality that allows velocity-dependent models to be mapped onto velocity-independent models in a well-defined way using the scattering time-scale computed when the halo achieves its minimum central density. We show how this time-scale depends on the halo parameters and an average cross-section computed at the central velocity dispersion when the central density is minimum. The predicted collapse time is fully defined by the scattering time-scale, with negligible variation due to the velocity dependence of the cross-section. We derive new self-similar solutions that provide an analytic understanding of the numerical results.